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1* INTRODUCTION 

We present here a steady-state, horizontally homogeneous, cloud-topped marine bound- 
ary layer model based primarily on the work of Lilly (1968) and Schubert et al. (1979). 
The conservative thermodynamic variables are equivalent potential temperature 9 e and to- 
tal water mixing ratio q + l. Some of the differences between this and Lilly’s (1968) model 
are: (1) radiation is allowed to penetrate into the boundary layer; (2) cloud top values 
of longwave radiation, equivalent potential temperature, and water vapor mixing ratio are 
linear functions of height derived from climatological data at California coastal stations; (3) 
the closure assumption assumes a weighted average of Lilly’s (1968) maximum and mini- 
mum entrainment theories. This model has been programmed in FORTRAN and will run 
interactively on an IBM-compatible personal computer. The program allows the user to 
specify the geographical location, the wind speed, the sea-surface temperature, the large 
scale horizontal divergence, and the initial guess for cloud top height. Output includes 
the steady state values of cloud top and cloud base height, mixed layer equivalent poten- 
tial temperature and total water mixing ratio, and the associated convective and radiative 
fluxes. The notation used throughout this abstract is that of Lilly with the exception that 
the subscript U replaces UH. 

2. RADIATION PARAMETERIZATION 

Since longwave cooling off cloud top is the important driving mechanism in the cloud- 
topped marine boundary layer, the radiation paramaterization warrants a detailed descrip- 
tion. We follow the philosophy that, consistent with the limitation of the vertical thermo- 
dynamic structure to two degrees of freedon, the vertical resolution of the radiative cooling 
should also be limited to two degrees of freedom. This means radiation can appear at most 
in the mixed layer thermal budget equation, and the cloud top jump condition. This is 
slightly more general than Lilly’s restriction of the radiative cooling to the cloud top jump 
condition. Thus, we write the equations for the changes in radiative fluxes across cloud top 
and across the mixed layer as 

Fu — Fh ~ ( P c p ) 1 1(1 - — c\j) — (1 — //)<!>| , (1) 

and 

Fh - Fs — {pc p )~ x {p(c iTfj - 4 ) - p's) , ( 2 ) 

where p is the constant air density, a is the Stefan-Boltzmann constant, T// is the cloud 
top temperature in Kelvin, 4 > s t he downward longwave radiative flux attributed to that 
portion of the atmosphere which lies above the mixed layer, S is the absorbed broadband 
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shortwave radiative flux, and p and p' are the longwave and shortwave radiation partitions 
respectively. These partitions can take on values from zero to one. If both are set to zero, 
the radiation parameterization reduces to Lilly’s case. Here and in the following sections 
the subscript U will refer to those properties of the upper air just above cloud top, the 
subscript H will refer to those properties just below cloud top, the subscript S will refer to 
those properties at the surface, and the subscript 0 will refer to those properties at some 
small height above the surface (typically taken to be 10 m). For this model, the downward 
longwave radiative flux term in (1) and (2) is expressed as a linear function of cloud top 
height. The function was derived based on the July average sounding data for the five years 
1976-1980 as extracted from the U.S. Dept, of Commerce Climatological Data, National 
Summary. Two California locations were used, San Diego and Oakland. The sounding data 
combined with the midlatitude ozone profile as taken from the U.S. Standard Atmosphere 
Table, 1976 and a uniform carbon dioxide profile of 0.501 gkg -1 was input into a broadband 
longwave radiation model written by Cox (1973). The model output was used to calculate 
the required linear relationships for the two locations. The resulting equations are given by 

c i _ f 314.0 + 0.03077/7 for Oakland, 

U 1 333.1 + 0.32360# for San Diego, ( 3 ) 


where the units for C\j and H are Wm 2 and m respectively. The value of the absorbed 
solar radiation is that suggested by Lilly (1968), i.e. S ~ 22.3 Wm -2 . 

To calculate X//, we start with the equation for potential temperature at cloud top (6h). 
The potential temperature at cloud top is equal to the potential temperature at cloud base 

(0 h ) plus the change that occurs when following a moist adiabat from cloud base to cloud 
top, i.e. 

>H = h + ; (iTz) (h - W 

In addition, 0 h can be calculated from our conservative thermodynamic variables by use of 
the definition of 6 e , i.e. 

9h = 9e - j-(q + l). (5) 

Cp 


To obtain the expression for the temperature at cloud top (T H ), we use both Poisson’s equa- 
tion and the integrated hydrostatic equation in potential temperature form. The resulting 
expression is 


where p H is determined hydrostatically from H. 


( 6 ) 


3. THE COMBINED CONVECTIVE-RADIATIVE MODEL 

With the discussion of the radiation parameterization completed, we can now write the 
combined convective-radiative model. The model equations, which are listed below, consist 
of two surface flux equations, the cloud base equation, the cloud top jump definitions, 
the cloud top temperature equation, the radiation equations, the consistency relation, the 
entrainment assumption, the mixed layer equivalent potential temperature and total water 
budget equations, and the cloud top jump condition on equivalent potential temperature. 


180 



These thirteen equations form a closed set in the thirteen unknowns if, ft, 6 e , (q + £), A 0 e , 
A (q + /), T„, ( Fu - F„), (F„ - Fs), (^)o, + 1') 0 , and w'(cf + £%. 


(V*e ) 0 = C T V(d eS - 0 e ) 


W(<t + l')o = C T V[qs-{q + e)] 

_ (1 + a)(qs ~ go) ~ a(0eS ~ #e) 

b 


h = 


A 6 e = 0 e (/ — 9 e 

A(q + £) = q v ~(q + £) 

Fu-F h = (pep)' 1 {(1 - p)(*T? { - C l v ) - (1 - //)<*} , 
Fh-F s = (pc p ) -1 {p(<rT^ - 4) - p's} , 


Th 


a \\ 

a 12 

.^21 

a 22 J . 


Me {W%h ~ ~ (F W ~ Fs) 

dt H 


djq + £) _ (wY) o - + 1% 

3/ H 


dH Fu-Fh- (w'0 e )„ 


dt 


Ad' 


DH. 


(7) 

( 8 ) 

(9) 

( 10 ) 

(ID 

( 12 ) 

(13) 

(14) 

(15,16) 

(17) 

(18) 
(19) 


Equations (15) and (16) are simply a shorthand matrix notation for the consistency relation 
and the closure assumption. In the following section we discuss the steady-state solutions 
to the above model. 


4. THE STEADY-STATE SOLUTIONS 


In the steady-state case all derivatives with respect to time are set to zero, which results 
in a system of nonlinear algebraic equations. The method we have chosen to use in solving 
for the steady-state solutions is to reduce our system of equations to one equation in H . 
We can then use a simple secant method algorithm to iteratively find its zero. Before 
desc ribing the steady-state model equations, it should be mentioned that the total water 
flux w , (q / -f £ f ) becomes constant with height. This can seen directly from (18). For this 
reason, the subscripts on this variable are deleted. 

The first step is to derive expressions for the surface fluxes of equivalent potential tem- 
perature and total water. The former is accomplished by eliminating the dependent vari- 
ables 0 e and (w'9 f e )H between the steady-state forms of (7), (17) and (19). The resulting 
expression is 


K*e)0 


(Fu - F h ) + (F h - F s ) + DH(0 eS - 6 eU ) 



DH \ 
c T v) 


( 20 ) 
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The latter expression is derived similarly, resulting in 


w'(q' + (') = 


DHjqs ~ qu) 


( 21 ) 


It should be noted that the cloud top values of equivalent potential temperature and water 
vapor mixing ratio (9 e jj and qu respectively) are linear functions of height derived from the 
same climatological data as used in the radiation paramaterization. 

We can also eliminate 9 e from our cloud base equation. This is accomplished by substi- 
tuting (7) and (8) into (9). The resulting expression is 


h _ (1 + a)w f (q / + £ f ) — fl(u/0p)o 

bC T V * 

To derive an expression for (w' W e ) H , we eliminate 0 e from (19) by using (7), which results 
in 

0?W e ) H = (F V - F„) - %)o + DH(8 eS - 9 eU ). (23) 

Finally, we need to derive the steady-state closure equation. Before this can be done, 
however, it is necessary to provide expressions for the virtual potential temperature flux at 
the surface, just below cloud base, just above cloud base, and at cloud top, respectively. 
These expressions expressions can be shown to take the form 


( W 'K ) o = K<?e)o - - 08 j w'(<f + l'), (24) 

( w ' e 'v)h - = {w%)h ~ (jr - ^ w'{</ + Fj, (25) 

(¥0i) h+ = 1+ 1 / 6 ) ) { ^ )h ~ + (26) 

= (~ + i+ 1 Q + ^ ) ) (*%)» - 9u/(g> + e% (27) 

where (w'9 f € )^ is given by 

(®^)a = (l - j) (^)o + (28) 


With the above variables defined, we can now derive our closure equation. We use a weighted 
average of Lilly’s (1968) maximum and minimum entrainment case to close our system 
equations. This closure equation takes the form 

k rH 

~H Jo W ' 6 'v dZ + ( J ~ = 0, (29) 

where A; is a weighting parameter which can take on values in the range 0 < k < 1. If we 
assume the daily averaged solar radiation is never strong enough to overcome the longwave 
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cooling and thus produce a net warming affect, the minimum virtual potential temperature 
flux must always be just below cloud base. This can be seen from (17). With this in mind, 
(29) can be integrated to give 

+ {w%)h + (w%)h+ + ^[(™%)s + ( w%)h - ~ (w r K)H ~ (w^C)fc+] = 0. 

( 3 °) 

We now have a closed set of equations (12)-(14), (20)-( 28), and (3 0) in the unknowns 
H, h, T h , ( F v - Fh ), (F h - F s ), (w’0' e ) H , { w'6 ' e ) h , w'{<f + t% (w'0' v ) H , (w'0' v ) h +, 

(«/<?' )/i-i and (w'0' v ) o . The method we have chosen to solve this system of equations can 
be explained in the following four step iteration sequence. 

(1) Make an initial guess of the cloud top height H and the cloud top jump in radiative 
flux (Fu - Fh). 

(2) Using the current estimates of H and (Ft/ — Fh), calculate in order (w'9' e ) o, w f (q f + V), 
h , (w'^h and the four virtual potential temperature fluxes using (20)-(27). 

(3) Again using the current estimate of 77, calculate the new radiation variables 7//, 
( Fjj — Fh ), and (Fh — Fs). This newly calculated value of (Fu — Fh) will be used in 
the next iteration if another iteration is required. 

(4) Using the above fluxes, check to see if (30) is satisfied to within some tolerable limit. 
If it is not, use the secant method to produce a new estimate of H and repeat steps 
(2)-(3) until (30) is satisfied. 

5. RUNNING THE PROGRAM 

The FORTRAN source code for the above model is stored on the attached floppy disk 
in the file STRATUS. FOR and the executeable code in the file STRATUS.EXE. To run the 
program, simply enter STRATUS. You will be prompted to enter the sea-surface tempera- 
ture, the large scale horizontal divergence, the wind velocity, and the initial guess for cloud 
top height. The output consists of cloud top and cloud base height, mixed layer equivalent 
potential temperature and total water mixing ratio, values of the convective and radiative 
fluxes and the cloud top temperature. 
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